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Abstract. We obtain factorized domain wall partition functions in two el- 
liptic height models: 1. A Felderhof-type model, which is new, and 2. A 
Perk-Schultz-type g£(l|l) model of Deguchi and Martin. 



0. Introduction 

0.1. Factorization in trigonometric vertex models. In [1], we obtained fac- 
torized domain wall partition functions (DWPF's j3 in two series of trigonometric 
vertex models: 1. The iV-state Deguchi- Akutsu models, for N G {2,3,4} (and 
conjectured the result for N > 5), and 2. The gl(r + 1 [s + 1) Perk-Schultz models, 
{r, s} G N (where given the symmetries of these models, the result is independent 
of r and s). 

0.2. Asymmetry. These models were characterized by an asymmetry of the vertex 
weights under conjugation of state variables. For example, in the Deguchi- Akutsu 
model, with state variables a £ {1, • • • , N}, the vertex weights are non-invariant 
under the conjugation a — > (N — a + 1). In the Perk-Schultz models, a similar 
property holds. Since one can trace the factorization of the DWPF's obtained 
in [1] to this asymmetry, it is natural to look for height modelfl with the same 
property. 

0.3. Factorization in elliptic height models. In this work we consider height 
models, where the state variables are heights, that live on the corners of the faces of a 
square lattice. A weight is assigned to each face. The weights are elliptic functions 
of the corresponding rapidities, external fields (if any) and height variables. As 
in [1], the models in this work are characterized by an asymmetry of the weights, in 
the sense that the weights of certain vertices (called line-permuting vertices) have 
different zeros, which leads to the factorization of the DWPF's. 

0.4. Summary of results. We obtain factorized DWPF's for two elliptic height 
models: 1. A Felderhof-type model, which (to the best of our knowledge) is new, 
and 2. A Perk-Schultz-type gl(l\l) model of Deguchi and Martin [3]. These are 
the first examples of DWPF's for elliptic and/or height models. 

0.5. Outline of paper. In Section 1, we collect a number of basic definitions 
related to elliptic height models to make the paper reasonably self-contained. In 
Section 2, we introduce a new Felderhof-type elliptic height model and obtain 
the corresponding factorized DWPF. In Section 3, we do the same for the <?/(l|l) 
elliptic height model of Deguchi and Martin. Section 4 contains brief remarks. 
The presentation is elementary in the hope that the paper will be reasonably self- 
contained. 
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1 For a review of previous results on the subject, see [2]. 

2 Also known as face or interaction-round-face (IRF) models. 
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1. Height models 

1.1. Faces and corners. We work on a square lattice, as in Figure [TJ with L 2 
square faces /y, where 1 < i < L increases from top to bottom, and 1 < j < L 
increases from left to right, /y has four corners that are labelled from top-left 
clockwise as {cjj, c it j + i,Ci + ij + i, (k+i,j}- 
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Figure 1. A square lattice with L 2 faces fij. Rapidities {u,v} and external fields 
{P; q} fl ow along lines that cross the faces. Height variables hij live on the corners. 



1.2. Heights and restrictions. We assign each corner cy a height variable hij, 
< i,j < L. In height models such as Baxter's solid-on-solid model [4], the heights 
are integral (possibly up to an overall shift). In the model of Section 2, the height 
variables depend linearly on the external fields which are continuous parameters, 
so they are no longer integral. We define the heights and the restrictions that they 
obey on a model by model basis in Sections 2 and 3. 

1.3. Flow lines, orientations and variables. There are L horizontal and L 
vertical lines that intersect at the middle points of fij. They indicate the flow of 
rapidities and external fields through fij. We assign the i-th horizontal line an 
orientation from left to right, a complex rapidity Ui and a complex external field 
Pi. We assign the j-th vertical line an orientation from bottom to top, a complex 
rapidity Vj and a complex external field qj, as in Figure [TJ 

1.4. Weights and Yang-Baxter equations. We assign each /y a weight Wy 
that depends on the height variables on its corners, the difference of the rapidity 
variables and the two external field variables (if any) flowing through it. The weights 
satisfy a set of Yang-Baxter equations. The weights and Yang-Baxter equations of 
the models discussed in this paper are given in Sections 2 and 3. 
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1.5. Elliptic functions and a theorem. Following the conventions used in [4], 
Chapter 15, we consider the elliptic function 

oo 

H{u) = J] (l-2q 2 «cos^+q 4 ") (l-q 2 ") (1) 

n— 1 

where u e C, q = cxp [ ^f- j , 21 and 21' (usually called 2K and 2K') are 

respectively the (real) width and height of an (upright) rectangle R in the complex 
u-plane, so that < q < 1. It is convenient to define 



1 J 2qV4 

which is entire and satisfies the quasi-periodicity properties 



(2) 



[u + 21} 



(3) 



r „ T /i 1 ( —Kin \ , . 
[ti + 2»J'] = --exp [ — J N 



(4) 



Theorem 1. If f(u) is an entire function that satisfies the quasi-periodicity con- 
ditions 



f(u + 2I) = (-) L f{u) 



(5) 



fiu + 2U>)={±) L ^[^pA) 



then 



/(«) 



(6) 



f L-\ 

/(«)=« nt"-^ 

where k and Ci, • • • , Cl-i are constants. 



L-l 

i=i 



(7) 



Proof. This is a refinement of Theorem 15(c) of [4], and the proof uses a similar 
argument. Choose the period rectangle R such that f(u) has no zeros on the 
boundary dR, and integrate 



/» 
/(«) 



on the anti-clockwise contour <9i?. From the 



quasi-periodicity conditions it follows that 



OH 



f'(u) 

/(«) 



(8) 



Hence the sum of residues of yj^y in R is equal to L, showing /(it) has exactly L 
zeros in i? (counting a zero of order n with multiplicity n). Writing the L zeros as 
Ci, • • • , (l, we define the function <j)(u) = Y\.j=i l u ~ G']- Since ^ log(/(u)/0(u)) is 
doubly periodic (by construction) and holomorphic (also by construction) one has 
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where A is a constant. Integrating, we obtain f(u) — ne Xu Y\j—\ [u — Cj]- Using 
the quasi-periodicity conditions of f(u), we can, without loss of generality, choose 
A = and Ql — V — 5Zj=i Oj which concludes the proof. 

2. A Felderhof-type height model 

In this section, we introduce an elliptic height model with weights that depend 
on rapidities, external fields and height variables. In the trigonometric limit, it 
reduces to the first in a series of models introduced by Deguchi and Akutsu in [5]. 
In that same limit, and decoupling the dependence on the height^, it reduces to 
the trigonometric limit of the elliptic Felderhof vertex model, which is the 2-state 
Deguchi- Akutsu model [6]. 

2.1. Notation. Given the rapidities {u, v} G C, external fields {p, q} G C, an 
upper-left corner height h e C, {Ai,A 2 } € {0,1}, A 3 G {0,1,2}, we use the 
notation 

Wuv f '; '; +q - A \ 1 (io) 

^ h+p-A 2 h+q+p-A 3 ) 

for the weight assigned to the verted represented in Figure [H 



h+q-A 1 



h+p-A 2 h+q+p-A 3 
I v,q 



Figure 2. A Felderhof-type face configuration. 



2.2. Height restrictions. For p = q = ^, we require that the heights satisfy the 
same restriction as in Baxter's solid-on-solid model [4], up to a normalization. More 
precisely, 



hij - h t+ i,j - ±-, hi : j - hij+i - ±- (11) 

2.3. The crossing parameter = |. The vertex weights will be parametrized in 
terms of the elliptic functions [it]. As defined in Equations [T1 and [2l [u] depends on 
the real parameters, I and which are the magnitudes of the half-periods of [u\. 
In the Felderhof-type model discussed in this section, we set / = 10. 



This can be achieved, for example, by introducing a parameter § G iR, shifting all height 
variables by § (the Yang-Baxter equations remain satisfied), then taking the limit f — > ioo. 
4 In the sequel, we simply say 'vertex' instead of 'face configuration'. 

^In the limit of zero external fields, that is p = q = i this is equivalent to setting the crossing 
parameter in Baxter's solid-on-solid model to the free fermion point. For details, see [4] 
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2.4. The weights. In the above notation, the non-zero weights are 



( h h+q \ 

Wuv { h+ P h+ g+P J 

w uv \ h 

y h+p-l h+q+p-2 ) 



a+(u,v,p,q) = [u-v+p + q] (12) 
a_(u, v,p, q) — [v - u + p + q] (13) 



W uv ( h h+p ^ ) = b + (u,v,p,q,h) (14) 

[2h]i[2{h + p + q)]i 



[2(h+p)]*[2(h + q)] 



t[u ~ v + q - p] 



W uv h = b-(u,v,p,q,h) 

^ h+p-l h+q+p-1 ) 

[2h]i[2{h + p + q)]i 



(15) 



[2(h + p)}i[2(h + q)} 



T [u-v+p-q] 



uu 



( h+q \ = c+(u,v,p,q,h) (16) 

h+p h+q+p-1 ) +V > y * 



[2(h + p)}i[2(h + q)\ 



T [v-u+p + q + 2h] 



W uv ( h h+p ^ ] = c.(u,v,p,q,h) (17) 

= r [u-v +P + q + 2h] 

[2(h+p)]h[2(h + q)}h l 

2.5. The Yang-Baxter equations. For rapidities {u, v, w}, external fields {p 7 q, r}, 
and non-negative integers {k, I, to, n, o}, the above weights satisfy the Yang-Baxter 
equations 



y w u A h h+q - ] 1 w uw \ h+q - j h+q+r - 1 1 

— ' I h+p— o h+q+p—n I 1 h+q+p— n h+q+r+p—m I 

{ h h+r-k \ 

x W vw 

\ h+q-j h+r+q-l J 

V w uv I h+r ~ k h+r+q - 1 }w u J h h+r - k } 

z — ' I h+r+p-j h+r+q+p-m I ; h+p-o h+r+p-j I 

j>Q v ' v ' 

x W vw f h+p -° h+p+r - ] 1 (18) 

I h+p+q— n h+p+r+q—m I 

Proof. This can be proved by direct computation using elliptic function identities, 
along the same lines as in [4]. For example, when {fc, Z, to, n, o} = {0, 1, 1, 1, 1}, the 
Yang-Baxter equation is 
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ft+ g+r+p— 1 



h 

h+p-1 



h+r 
h+r+p 



(19) 



Using the expressions for the weights, we obtain 



c-(u,v,p,q,h) a+(u,w,p,r) b-(v,w,q,r,h) 
+ b-(u,v,p,q,h) c-(u,w,p,r, h + q) c+(v,w,q,r,h) 
= c-(u,v,p,q,h + r) b-(u,w,p,r,h) a+(v,w,q,r) (20) 

Writing the weights in terms of elliptic functions, one can eliminate common 
factors, and the proof of the equation reduces to proving 



[it — v + p + q + 2h] [it — w + p + r][v — w + q — r] [2(h + q + r)] 
+ [u - v + p - q] [u - w + p + r + 2{h + q)] [w - v + q + r + 2h] [2r] 
= [u~v+p + q + 2(h + r)}[u-w+p-r}[v-w + q + r}[2(h + q)} (21) 

which proceeds by noting that the ratio of the left-hand-side and right-hand-side is 
doubly periodic and entire in u, and therefore a constant with respect to u. Setting 
it = v — p + q, the constant is found to be 1. 



2.6. Switching off the external fields. Setting p = q = | is equivalent to 
switching off the external fields. This becomes clear by inspection of the vertex 
weights, which up to normalization become equal to those of Baxter's solid-on- 
solid model at the free fermion point. 



2.7. The external fields tilt the heights. One can think of the external fields 
p ^ \ and/or q ^ |, as effectively tilting the heights of the lattice faces that they 
flow through. This tilt is with respect to the line along which a field flows. This 
effectively adds to or subtracts from the height differences that are the case in the 
absence of external fields. 



2.8. The c + vertex. In discussions of DWBC's and DWPF's, the c+ vertex, see 
Figure [3l plays a special role: It is the DWPF for a 1 x 1 square lattice. 
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h+p 



h+q+p-1 



I «,5 



Figure 3. TTie Felderhof-type c + vertex. 



2.9. Domain wall boundary conditions (DWBC). We define the DWBC's as 
an expanded c+ vertex, as in Figure [4] Given the external fields {p, q} and starting 
from /ioo — h at the top-left corner, the boundary heights change by qj from left to 
right along the upper boundary, pi — 1 from top to bottom along the right boundary, 
— Qj + 1 from right to left along the lower boundary, and — pi from bottom to top 
along the left boundary. 
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h + qi.L +Px-1 



h + qi.L +Pi,2 - 2 



h + qi.L +Pi.l - L 



Figure 4. Felderhof-type height domain wall boundary conditions. We use the 
notation p^ 3 = Y^i=iPk, etc - 



2.10. Domain wall partition function (DWPF). The DWPF on an L X L 
lattice, Zlxl, is the sum over all weighted configurations that satisfy the DWBC. 
The weight of each configuration is the product of the weights of the vertices 
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z L*.L — 

configurations 



n 



(22) 



2.11. Line permuting vertices. In proofs of DWPF's two vertices play an im- 
portant role. These are the o-type vertices which can be used to permute adjacent 
flow lines. 



u, p 



h+p 



h+q 



h + q-1 



u, p 



h+q+p 



h+p-1 



h+q+p-2 



T 

I v,q 



T 

I v,q 



Figure 5. The Felderhof-type line permuting vertices a + and a_. 

2.12. Different zeros. The weights of the line permuting vertices, [u — v + p + q] 
and [v — u +p + q], have different zeros. This is the property that will allow us to 
obtain the zeros of the DWPF and compute it in factorized form. 

2.13. Properties of the partition function. The following four properties de- 
termine the partition function uniquely. 

2.13.1. Property 1: Quasi-periodicity. The partition function is entire in ui and 
satisfies the quasi-periodicity conditions 

Zlxl [u 1 + 2,...,u L ,{v},{p},{q},h] =(-) L Z LxL [{v},{v},{p},{q},hj 

(23) 



Zlxl - 2l l °^ q \ {v}, M, {<?}, h j 



— iri 



Lu x + {L- 2) Pl - ^2( Vj +qj)- 2h 



Zl^l [{u},{v},{p},{q},h} (24) 

Proof. Since the weights are entire functions in the rapidities, it follows that 
Zlxl ^ {u}, {v}, {p}, {q}, h j is an entire function in u\. To prove the quasi- 
periodicity conditions, we write the partition function in the form 



Zlxl [{u},{v},{p},{q},h} = ^ P n [u u {v}, Pl ,{q},h] x (25) 

n=l 

Q n \u 2 , . . . , u L , {v},p 2 , ...,p L , {q}, 
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where 



{v},Pi,{q},h} = 



k j=n+l fe=l 



l,«n,Plign,fe + y^gfc 



fc=l 



and Q„ ^«2, • • • , ul, {v},P2, ■ ■ ■ ,Pl, {(?}, h j does not depend on u\. Using the 
expressions for the weights, we have 

P n (ui+ 2, {«},pi, ) ={-) L Pn («l,W,Pl, {«},/») 



■ exp 



-TTl 



(27) 
(28) 



Ltii + (L - 2)p! - ^2(vj +qj)- 2h 

3 = 1 



Pn [u 1 ,{v},p 1 ,{q},h} 



from which the required property follows immediately. 



2.13.2. Property 2: Simple zeros. The partition function has simple zeros at u\ = 
uj — pi — pj, where j = 2, . . . , L. 

Proof. We multiply the partition function by a + (u2, U\,p2,pi), and use the Yang- 
Baxter equation to slide the inserted face through the lattice. 



U 2 ,P2 



Ml, Pi 




Ml, Pi 



U2,P2 



Figure 6. Inserting an a + vertex into the left boundary. 



It emerges as a_(it2, Ui,p2,Pi), and the order of the first two lattice rows is 
reversed. 
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"2,f>2 



Ul,Pl 




Ml, Pi 



U 2 ,P2 



Figure 7. Extracting an a_ vertex from the right boundary. 



This is equivalent to the equation 



Z L y.L [{u},{v},{p},{q},h} = 



a-(u 2 ,ui,p2,pi) 



(29) 



a + (u 2 ,u 1 ,p 2 ,p 1 ) 

ZlxL \u2,Ui,.. .,U L ,{v},p 2 ,Pl, --^PL, {q},h^ 

Repeating this procedure on the second and third rows, and so on, we obtain 

Zl*l ({m},W,W,M^) = f[ \ a -^>*>*\ } x (30) 

Z L y.L ["2, • • .,u L ,u 1 ,{v},p 2 , . . . ,Pl,Pi, {q}, h j 
which has the required simple zeros in the numerator. 

2.13.3. Property 3: A recursion relation. The partition function satisfies the recur- 
sion relation 

Z L y.L \{u},{v},{p},{q},h\\ =c+(v 1 -p 1 -q 1 ,v 1 ,p 1 ,q 1 ,h)x 

\ J \u\=v\—p\-q\ 



( L 



Y[b + (u j ,v 1 ,p j ,q 1 ,h + '^2pk) b-(vi - pi - qi, v j ,p 1 , q jt h + q k ) 



fc=l 



fe=l 



Z(l-i)x(l-i) ("2, ...,u L ,v 2 ,.. -,v L ,p2, ■ ■ ■ ,PL,q 2 , ...,qL,h + pi + qi-l \ (31) 

Proof. In any lattice configuration in the partition function sum, the top-left corner 
of the lattice must be a + (u\, V\,pi, qi) or c + (u\,vi,pi,qi,h). Setting u\ — v\ —p\ — 
qi in the partition function sets to zero all configurations with a+(ui, Vi,pi, q\). 

The surviving configurations must have a top- left corner equal to c + , which fixes 
the rest of the top row to 6_ , the rest of the first column to b+ , and the remainder 
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of the lattice to Z(l-i)x(l-i)- The above recursion follows immediately from these 
considerations. 

2.13.4. Property 4- The partition function on a 1 x 1 lattice is given by 

Z lxl (u 1 ,v 1 ,p 1 ,q 1 ,h) = c + (u 1 ,v l7 p 1 ,q 1 ,h) (32) 
Proof. This follows from the definition of domain wall boundary conditions. 

2.14. The partition function is uniquely determined. Assume that Zt n _i\ x t n _u 

is uniquely determined by the above four properties, for some n > 2. From Property 
1, Property 2 and Theorem 1.1, we have 



({«}. M,{P}> {?}.>*) = K(u 2 ,...,u n ,{v},{p},{q}) 



E( 



)+S(pj+?j)+2/»: 



JJ[ttl -Mj +P1 

U=2 



(33) 



Property 3 fully determines the coefficient re in terms of ■^( n -i)x(n-i)- Finally, 
since Zi x i is uniquely determined by Property 4, Z nxrt is uniquely determined by 
the four properties. 



2.15. The domain wall partition function. The solution to the preceding four 
properties is given by 



v il \ S \ S \ S \ h\ Ili=i[2Pj]*[2<&]* 
Zlxl [{u},{v},{p},{q},hj 



[2(^ + E^i^)] 5 [2(^ + Et=i*)] 3 

L L 

E ~ u 3 ) + ^ + * ) + 2/i] ]J t M J ~ Mfc + p J + t Wfe ~ ^' + ?fe + 9j] 

j=l j = l l<j<k<L 

(34) 



2.16. Two choices of DWBC's. There are two possible choices of DWBC's. One 
corresponds to an expanded c+, as in this work, and one to an 'expanded' c_ vertex. 
The DWPF depends on the choice. The two expressions coincide for vanishing 
external fields, that is pi = qj = ^, and appropriate choices of the boundary height 
variables. 



3. A Perk-Schultz-type g/(l|l) height model 

In [3], Deguchi and Martin introduced elliptic height versions of the gl(r+l\s + l) 
trigonometric vertex models. In the following, we define DWBC's and compute the 
DWPF in the gZ(l|l) case. Since the analysis in this Section follows almost verbatim 
that of Section 2, we will be brief and give just enough details where the two models 
differ. 
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3.1. Notation, heights and restrictions. In this model, there are two square 
lattices. A physical L x L lattice that the heights live on, and a target 1x1 lattice 
that the heights take values in. The target lattice is spanned by the unit vectors e M , 
fi 6 {— 1, +1}, thus the height variables are 2-component vectors /i„}. Heights 
on adjacent physical lattice corners are restricted to take values in adjacent points 
on the target lattice. Height differences along e_i and e+i lead to different vertex 
weights. To each height vector h = {h^, we assign a scalar 



(35) 



where u> is an arbitrary constant antisymmetric complex 2x2 matrix. We use 
as well as e s i gn ^ to indicate the same unit vector. 

3.2. No external fields and the crossing parameter is a variable. Unlike 
the previous Felderhof-type model, the Perk-Schultz-type model in this Section has 
no external fields. The crossing parameter is left as a variable. 



3.3. The weights. The non-zero vertex weights are 



W uv 



h 

h + e^ 
h 

h + e. 



h 



h + 2e, 



= a M (u,w) = 



[1 + fi(u — v) 



(36) 



h + e M + e h 



= b f _ l (u,v) = MUL 1 , fiftv 



[hp V - (it - v)} 



(37) 



V.±v (38) 



The variables huv on the right hand sides of Equations 1371 and 1381 are the scalars 
assigned to the heights at the upper left corners of the corresponding vertices. 

3.4. The Yang-Baxter equations. The weights satisfy the following Yang-Baxter 
equations [3]. 



) W ™ [Is] W ^ ( 



c 
d 

d 

e 



(39) 



3.5. The c + vertex. We take the c + vertex to be as in Figure [8j 



I1+C4 



h+ex. -f-e- 



Figure 8. A general face configuration. 
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h + e+ 
h + 2e 4 



h + Le A 



h + Le + 



h + Le + + e_ 
h + Le+ + 2e_ 



h + Le_ + Le.A 



f t 

h + Le + + 2e_ 
h + Le + + e_ 

Figure 9. Perk-Schultz-type height domain wall boundary conditions. 



3.6. The domain wall boundary conditions. We choose the DWBC as in Fig- 
ure [9] In other words, starting from the lower left corner, all height changes along 
the left boundary are of type e+i, along the upper boundary they are of type e_i, 
along the right boundary they are of type e+i, then along the lower boundary they 
are of type e-%. 

3.7. The line-permuting vertices. We take the line permuting vertices to be as 
in Figure 1101 Their weights have different zeros leading to a factorization of the 
DWPF, just as in Section 2. 



h+e + 



h+e+ 



h+e. 



h+2e+ 



h+e- 



h+2e- 



i, 



Figure 10. The Perk-Schultz-type line permuting vertices a+ and a_ . 



3.8. The DWPF. Having defined the model, the derivation of the corresponding 
DWPF proceeds precisely in analogy with that in Section 2. Based on the quasi- 
periodicity properties of the partition function, we propose a factorization in terms 
of the [u] functions. We obtain the zeros by permuting adjacent flow lines, using 
the line permuting vertices, obtain a recursion relation the DWPF satisfies and an 
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initial condition. The DWPF is uniquely determined, and the following expression 
satisfies all the conditions. 



[h + (L - 1) - Y,k=i( u k ~ Vk ^ TT [1 + Ui - Uj] [1 - (vi ~ Vj)] 
[h+(L-l)] 11 s [i] [i] 

(40) 



3.9. On the gl(r + l\s + 1) height models. The reason we chose the <?Z(1|1) 
case is that, given the way that we define DWBC's, the symmetries of the gl(r + 
l\s + 1) models are such that only the two state variables variables that we put 
on the domain wall boundaries end up propagating inside the configurations. This 
effectively restricts the DWPF to the gl(l\l) model. 

4. Remarks 

The point of this work is to give examples of domain wall partition functions in 
elliptic and/or height models. We restricted our attention to models that are non- 
invariant under state variable conjugation, which greatly simplified the problem. 

Because of the fermionic nature of the models discussed in this paper, there is 
no interesting applications of their DWPF's to enumerations of alternating sign 
matrices or related objects that we are aware of. 

It is highly likely that all models that are non-invariant under some form of state 
variable conjugation, that allows factorization, are fermionic (as the Felderhof-type 
model of Section 2) or contain fermions that play an essential role in the definition of 
the DWBC's (as the Perk-Schultz-type model of Section 3). As such, these models 
are unrepresentative of the general case. However, they are non-trivial and we 
hope that one can learn something by extending our results to compute correlation 
functions. 
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